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Abstract 

We report the measurement of a new observable of jet quenching in central Pb-Pb collisions at ^snn 
= 2.76 TeV, based on the semi-inclusive rate of charged jets recoiling from a high transverse mo¬ 
mentum (high-pj) charged hadron trigger. Jets are measured using collinear-safe jet reconstruction 
with infrared cutoff for jet constituents of 0.15 GeV/c, for jet resolution parameters R = 0.2, 0.4 and 
0.5. Underlying event background is corrected at the event-ensemble level, without imposing bias on 
the jet population. Recoil jet spectra are reported in the range 20 < Pt jet < lOOGeV/c. Reference 
distributions for pp collisions at yC = 2.76 TeV are calculated using Monte Carlo and NLO pQCD 
methods, which are validated by comparing with measurements in pp collisions at i/s = 7 TeV. The 
recoil jet yield in central Pb-Pb collisions is found to be suppressed relative to that in pp collisions. 
No significant medium-induced broadening of the intra-jet energy profile is observed within 0.5 radi¬ 
ans relative to the recoil jet axis. The angular distribution of the recoil jet yield relative to the trigger 
axis is found to be similar in central Pb-Pb and pp collisions, with no significant medium-induced 
acoplanarity observed. Large-angle jet deflection, which may provide a direct probe of the nature of 
the quasi-particles in hot QCD matter, is explored. 
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1 Introduction 

Hadronic jets are unique probes of the hot Quantum Chromodynamic (QCD) matter generated in nuclear 
collisions at collider energies. Interactions of hard-scattered partons with colored matter may modify 
intra-jet structure, softening and broadening the distribution of hadronic jet fragments relative to jets 
generated in vacuum, and may deflect jets by large angles. These phenomena, known as jet quenching 
ill , can probe dynamical properties of the hot QCD medium lO and the nature of quasi-particles in the 
Quark-Gluon Plasma (QGP) @. 

Jet quenching generates marked, experimentally observable effects. Measurements of inclusive distribu¬ 
tions and correlations of high transverse momentum (high-px) hadrons have revealed significant yield 
suppression in nuclear collisions relative to vacuum lf4l- IT6ll . Suppression of the inclusive yield of recon¬ 
structed jets Ifl7lj20ll and enhancement in the rate of energy-imbalanced back-to-back di-jet pairs B2lll22t 
have also been observed in nuclear collisions. A measurement of event-averaged missing p-\ suggests 
that the radiated energy induced by the interaction of an energetic parton with the medium is carried to a 
significant extent by soft particles at large angles relative to the jet axis f23l . 

The measurement of reconstructed jets over a wide range in jet energy and jet resolution parameter (R) is 
required for comprehensive understanding of jet quenching in heavy-ion collisions. Such measurements 
are challenging, however, due to the presence of complex, uncorrelated background to the jet signal, 
and the need to minimize biases in the selected jet population imposed by background suppression tech¬ 
niques. Multiple, complementary measurement approaches, differing both in instrumentation and in 
analysis algorithm, are therefore important to elucidate the physics of jet quenching using reconstructed 
jets. 

In this article we present a new approach to the measurement of jet quenching, based on the semi- 
inclusive distribution of charged jets recoiling from a high-px charged hadron trigger (“h-jet” coinci¬ 
dence) in central (0-10%) Pb-Pb collisions at ^/.vnn = 2.76 TeV. Jets are reconstructed using charged 
particle tracks with the kj Il24l and anti-kx algorithms G51 . with infrared cutoff for tracks px,const >0.15 
GeV/c. Uncorrelated background to the recoil jet signal is corrected solely at the level of ensemble- 
averaged distributions, without event-by-event discrimination of jet signal from background, using a 
technique that exploits the phenomenology of jet production in QCD. The correction is earned out using 
an unfolding technique. This approach enables the collinear-safe measurement in heavy-ion collisions of 
reconstructed jets with low infrared cutoff over a wide range of jet energy and R. Recoil jet distributions, 
which are differential in ;?x,jet and in azimuthal angle relative to the trigger axis, are reported for R = 0.2, 
0.4 and 0.5, over the range 20 < Pxjet < 100 GeV/c. 

Suppression of the recoil jet yield due to quenching is measured by comparison to the yield in pp col¬ 
lisions. However, our current data for pp collisions at = 2.76 TeV do not have sufficient statistical 
precision to provide a reference for the Pb-Pb measurements reported here. The reference distribution is 
therefore calculated using the PYTHIA event generator l26l and perturbative QCD (pQCD) calculations 
at Next-to-Leading Order (NLO) I l27l l. which are validated by comparison with ALICE measurements 
of pp collisions at sfs = 7 TeV. Angular broadening of the internal jet structure due to quenching is in¬ 
vestigated by comparing the differential recoil jet distributions for different values of R. Acoplanarity 
between the trigger hadron and recoil jet directions is measured to explore the deflection of the jet axis 
induced by quenching. The rate of large angular deviations is measured; this rate may be dominated by 
single hard (Moliere) scattering, which could potentially probe the quasi-particle nature of the hot QCD 
medium 011281. 

These observables are directly comparable to theoretical calculations, without the need to model the 
heavy-ion collision background, due to utilization of a hadron trigger, the semi-inclusive nature of the 
observables, and the background suppression technique. The only non-perturbative component required 
to calculate the hard-process bias is the inclusive charged hadron fragmentation function (in-vacuum or 
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quenched) for the trigger hadron. 

The paper is organized as follows: Sect. [2j dataset, event selection, and simulations; Sect. [3j jet recon¬ 
struction; Sect. [TJ discussion of observables; Sect. 0 raw distributions; Sect. [6j corrections; Sect. [7] 
systematic uncertainties; Sect. [8j reference distributions for pp collisions; Sect. [9j results; and Sect. flOl 
summary. 

2 Data set, offline event selection, and simulations 

The ALICE detector and its performance are described in 12911301 . 

The Pb-Pb collision data were recorded during the 2011 LHC Pb-Pb run at ^/snn = 2.76 TeV. This 
analysis uses the 0- 10% most-central Pb-Pb collisions selected by the online trigger based on the hit 
multiplicity measured in the forward VO detectors. The online trigger had 100% efficiency for the 0-7% 
interval in centrality percentile, and 80% efficiency for the 8 - 10% interval. 

Events are reconstructed offline as described in Ref. flTl . Charged tracks are measured in the ALICE 
central barrel, with acceptance |rj| < 0.9 over the full azimuth. Accepted tracks are required to have 
0.15 < pj < 100 GeV/c, with at least 70 Time Projection Chamber (TPC) space-points and at least 
80% of the geometrically findable space-points in the TPC. To account for the azimuthally non-uniform 
response of the Inner Tracking System (ITS) in this dataset, two exclusive classes of tracks are used lf30ll : 
tracks with Silicon Pixel Detector (SPD) hits (70% of all tracks in central Pb-Pb collisions, and 95% in 
pp collisions); and tracks without SPD hits but with a primary vertex constraint. The primary vertex is 
required to lie within 10 cm of the nominal center of the detector along the beam axis, and within 1 cm 
of it in the transverse plane. After offline event selection, the Pb-Pb dataset consists of 17M events in 
the 0- 10% centrality percentile interval. 

The pp collision data used to validate PYTHIA and pQCD calculations were recorded during the 2010 
low-luminosity pp run at ^/s = 7 TeV, using a MB trigger. The MB trigger configuration, offline event 
selection, and tracking are the same as described in OTll . After event selection cuts, the pp dataset consists 
of 168M events. There is negligible difference in the inclusive jet cross section for events selected by the 
ALICE online trigger, and for a non-single diffractive event population. 

Simulations of pp collisions were carried out using PYTHIA 6.425, with the Perugia 0, Perugia 2010, and 
Perugia 2011 tunes Il32l . Instrumental effects are calculated using the Perugia 0 and Perugia 2010 tunes 
for pp and Pb-Pb collisions respectively, with a detailed detector model implemented using GEANT3 
Il33l . In addition, a simulation based on HIJING If34l is used to evaluate the detector response in the high 
multiplicity environment of Pb-Pb collisions. Perugia 2011, which has been tuned to other LHC data, 
is used as an alternative to compare with the new data presented here. Simulated events, which include 
primary particles and the daughters of strong and electromagnetic decays but not instrumental effects or 
the daughters of weak decays, are denoted “particle level”. Simulated events also including instrumental 
effects and weak decay daughters where reconstructed tracks are selected using the experimental cuts are 
denoted “detector level”. 

For central Pb-Pb collisions, tracking efficiency is 80% for pj > 1 GeV/c, decreasing to 56% at 0.15 
GeV /c. Track momentum resolution is I % at pj = I GeV/c and 3% at pj = 50 GeV/c. For pp collisions, 
tracking efficiency is 2%-3% higher than in central Pb-Pb collisions. Track momentum resolution is 1% 
at pj = I GeV/c for all reconstructed tracks; 4% at pj = 40 GeV/c for tracks with SPD hits; and 7% at 
p T = 40 GeV/c for tracks without SPD hits 13011311 . 
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3 Jet reconstruction 

Jet reconstruction for both the pp and Pb-Pb analyses is carried out using the kj Il24l and anti-kx 1251 
algorithms applied to all accepted charged tracks. The boost-invariant pj -recombination scheme is used 
If24l . Jet area is calculated by the Fastjet algorithm using ghost particles with area 0.005 1 1351 . 

Charged jets are not safe in perturbation theory, because radiation carried by neutral particles is not 
included. However, infrared-safe calculations of charged-jet observables can be performed using non- 
perturbative track functions, which absorb infrared divergences and describe the energy fraction of a 
parton carried by charged tracks l l36l . Track functions are analogous to fragmentation functions, with 
DGLAP-like evolution, and perturbative calculations using them are in good agreement with PYTHIA 
calculations |[36Tl . Track functions can provide the basis for rigorous comparison of the charged-jet 
measurements reported here with both analytic and Monte Carlo QCD calculations. 

For the Pb-Pb analysis, adjustment of jet energy for the presence of large background utilizes the FastJet 
procedure 071 . in which jet reconstruction is carried out twice for each event. The first pass applies the 
kj algorithm with R = 0.4 to estimate p, the density of jet-like transverse-momentum due to background, 
which is defined as 


{ raw,i A 


where index i runs over all jet candidates in an event, and an< ^ ^j e t are transverse momentum 
and area of the i th reconstructed jet. Further details are presented in B38l . The central data points in this 
analysis are determined by excluding the two jets with highest jet fr° m calculation of the median, with 
a variant used to study systematic sensitivity to this choice. 

The second pass, which generates jet candidates for the reported distributions, applies the anti Ay algo¬ 
rithm with resolution parameter R = 0.2, 0.4, and 0.5. The value of p^jet ' s correcte d according to 11371 . 


reco,i raw.i _ A i 
P T ,jet =PTjet -P' A V 


jet’ 


( 2 ) 


where and A- et are measured for the i th jet in an event, and p is a scalar value common to all jets in 
each event, but varies from event to event. 

Jet candidates are accepted if |rjj e t| < 0.5 for R = 0.2 and 0.4 and 177j et | < 0.4 for R = 0.5, where r/j et is the 
pseudo-rapidity of the jet candidate centroid. The azimuthal acceptance of the recoil yield measurement 
is n — A (p < 0.6, where A (p = <p ln „ — <pj et | is the difference between the azimuthal angle of the trigger 
hadron (<p tr ig) and the jet candidate centroid (<pj et ), and 0 < A<p < n. 

A cut on jet area is applied to suppress combinatorial jets while preserving high efficiency for true hard 
jets (39]|40|]. Jet candidates are rejected if 4j et < 0.07 for R = 0.2; 4j et < 0.4 for R = 0.4; and A- et < 0.6 
for R = 0.5. 

Similar procedures are followed for the pp data analysis. Jets are reconstructed with the anti-kj algorithm 
for R = 0.2, 0.4 and 0.5. Reconstructed Py j° t ' ch is adjusted using Eq. [2] where p is estimated in this case 
by the summed p-\ in two cones of radius R = 0.4, with centroids at the same i] but perpendicular in 
azimuth to the leading jet in the event. 

The instrumental jet energy resolution (JER), which characterizes the detector response relative to charged 
jets at the particle level, varies from 20% at PT,jet = 20 GeV/c to 25% at PT,jet = 100 GeV/c, for both 
Pb-Pb and pp collisions, with negligible dependence onR. The jet energy scale (JES) uncertainty, which 
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is dominated by the uncertainty of tracking efficiency, is approximately 5% for both Pb-Pb and pp col¬ 
lisions, with negligible dependence on pj h j eL and R. However, the instrumental response is significantly 
non-Gaussian HU and unfolding of the full response matrix is used for corrections. 


4 Discussion of observables 

4.1 General considerations 

Energetic jets that arise from high momentum transfer (high-<2 2 ) scattering of partons are readily visible 
in event displays of high multiplicity heavy-ion collisions 112211231 . However, accurate measurement of 
jet energy in such events, and unbiased measurement of jet distributions, are more difficult. Application 
of a jet reconstruction algorithm to high multiplicity events will cluster hadrons arising from multiple 
incoherent sources into each reconstructed jet, resulting in significant smearing of the true hard jet energy 
distribution. It will also generate a large population of “combinatorial” background jets comprising solely 
hadrons generated by soft production processes ( Q 2 below a few GeV 2 ), which cannot be identified as 
hard jets with smeared energy. 

Current heavy-ion jet analyses select the hard jet population on a jet-by-jet basis by several different 
approaches: removal of an estimated background component of transverse energy prior to jet reconstruc¬ 
tion BTil : or imposition of a fragmentation bias requiring a cluster of high-px tracks or a single high-px 
track in the jet, and imposition of a jet px threshold IU71ll9il20ll42ll . These rejection techniques may bias 
towards certain fragmentation patterns in the accepted hard jet population. 

This analysis takes a different approach, in which corrections for background and instrumental effects 
are applied solely at the level of ensemble-averaged distributions, without rejection of individual jet 
candidates or removal of event components. The analysis is based on the semi-inclusive differential 
distribution of charged jets recoiling from a high-px trigger hadron, with the trigger hadron selected 
within a limited px.trfg interval (Trigger Track, or TT, class). This distribution, which is the number of 
jets measured in the recoil acceptance normalized by the number of trigger hadrons, is equivalent to the 
ratio of inclusive production cross sections, 


d2N AA 


dp^ jet dT7 jet 


1 


£j2^yAA—^h+jet+X \ 


7 AA^h+X 


PT.trigeTT 


dp! 


Tjetdfijet 




PT.hGTT 


(3) 


where AA denotes pp or Pb-Pb collisions, ct aa >h ^ x is the cross section to generate a hadron within 
the px interval of the selected TT class, d 2 a AA ^ h+ j et+x /dpj!j et dr] is the differential cross section for 
coincidence production of a hadron in the TT interval and a recoil jet, and py h j et an< 3 r/jet are the charged 
jet transverse momentum and pseudo-rapidity. 

Because the observable in Eq. [3] is semi-inclusive, the selection of events containing a hard process 
(“hard-process selection”) is based solely on the presence of a high-px hadron trigger. In particular, 
there is no requirement that a jet satisfying certain criteria be found in the recoil acceptance. Rather, 
all jet candidates in the recoil acceptance are counted in Eq. [3} regardless of their specific properties. 
Events with no hard jet candidates (however defined) falling within the acceptance are not rejected, and 
contribute to the normalization. This observable thereby measures the absolutely normalized rate of 
recoil jets observed per trigger. Correction for the contribution of uncorrelated background jets in Eq. [3] 
is earned out at a later step in the analysis, as discussed below. 

Other jet correlation measurements in heavy-ion collisions have been carried out, in which hard-process 
selection utilizes a compound condition that requires the presence of both a trigger object (jet or photon) 
and a recoil jet satisfying certain criteria I22ll23ll43l . The jet correlation distributions in these analyses 
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are normalized per trigger-recoil pair; absolute normalization requires scaling by the inclusive trigger 
yield, together with selection of the recoil jet population using the semi-inclusive procedure described 
above. The role of normalization in the measurement of in-medium large-angle scattering is discussed 
in Sect. 19.21 

4.2 Trigger hadrons and hard-process bias 

The use of high-px hadron triggers for hard-process selection in this analysis is based on the following 
considerations. 

Hadrons with p-\ larger than about 5-7 GeV/c are expected to originate primarily from fragmentation 
of energetic jets, in both pp and Pb-Pb collisions at ^.vnn = 2.76 TeV (see e.g. |[44ll ). They provide 
experimentally clean triggers, without the need for correction for uncorrelated background. Selection of 
events by requiring the presence of a high-px hadron biases towards events containing a high-Q 2 partonic 
interaction, with jets in the final state. 

Inclusive distributions of high-px hadrons have been measured and calculated theoretically in both pp 
and heavy-ion collisions at collider energies. For pp collisions at the LHC, agreement within a factor 
two is found between NLO calculations and data for px > 10 GeV/c, with the discrepancies attributable 
to poorly constrained gluon fragmentation functions that can be improved by fitting to LHC data 1451 . 
For inclusive hadron production in heavy-ion collisions ll5~i?9l! 1 TI1T4TI . the medium-induced modification 
and evolution of fragmentation functions have been calculated in several frameworks, showing good 
agreement with data ( (4611471 and references therein). Suppression of inclusive hadron production in 
heavy-ion collisions has been used to determine the jet transport parameter q 0. 

Any hard-process selection procedure imposes bias on the accepted event population, and accurate com¬ 
parison of theory calculations with such measurements requires calculation of this selection bias. In this 
analysis, hard-process selection uses the same cuts as those used for high-px inclusive hadron measure¬ 
ments. Since inclusive hadron production is calculable in both pp and Pb-Pb collisions, the selection 
bias in this analysis is likewise calculable using current theoretical approaches. 

4.3 Hadron-jet coincidences 

There are additional considerations for interpreting hadron-triggered recoil distributions in Eq. [3] and 
comparing such measurements to theoretical calculations, as follows. 

The h-jet coincidence cross section in pp collisions has been calculated in a pQCD framework (27l . In 
this process at FO, a pair of final-state partons is generated with opposing transverse momenta, with 
one of the pair fragmenting into a hadron which carries momentum fraction z of the recoiling jet. Since 
z = px,trig/ pt, jet < 1 at FO, the requirement of a high-px,trig hadron above threshold therefore biases 
against coincident recoil jets with px.jet < Px.trig> but does not impose a kinematic constraint on recoil 
jets With Px.jet > Pt. trig ■ 

The inclusive hadron distribution at high-px is biased towards high-z jet fragments, due to inteiplay 
between the shape of the inclusive jet px spectrum and the shape of the inclusive fragmentation function, 
with (zind) ~ 0.6 at LHC energies 11481 . However, in the semi-inclusive measurement based on Eq. [3] the 
trigger-normalized rate of recoil jets is measured as a function of ■ At LO this corresponds to z = 
/A,trig/Px,jet> which can differ significantly from (zinci) (HI. For fixed px.uig the z-bias therefore varies 
as a function of recoil px.jet, with stronger bias than the inclusive case for px.jet ~ px.n-ig, and weaker 
bias for px.jet Px.trig- The z-bias has been calculated for pp collisions at y/s = 2.76 and 7 TeV using 
the approach of 1271 and found to be similar at FO and NFO. This z-bias, which is kinematic in origin, 
likewise occurs in nuclear collisions in which jets experience quenching. This effect is intrinsic to any 
theoretical framework based on pQCD, both in-vacuum and quenched, and will be properly accounted 
for in such calculations. 
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For quenched jets in nuclear collisions, high-/ j-j hadron selection may generate two additional, related 
biases. The first is a bias towards high-z fragments of jets that have lost relatively little energy in the 
medium |[49l . The second is a geometric bias in which small energy loss corresponds to small path 
length in matter. In the latter case, jets generating high-px trigger hadrons are generated predominantly 
on the surface of the collision region and headed outward OOl456ll . with the corresponding recoil jet 
population biased towards longer path length in matter than the unbiased, fully inclusive jet population. 

The degree to which high-px hadron selection biases towards small energy loss of its parent jet de¬ 
termines the degree of similarity in the underlying distribution of high-<2 2 processes in pp and Pb-Pb 
collisions, for the same hadron trigger cuts. In the following, we refer to potential differences in such 
Q 2 distributions as being due to “trigger-jet” energy loss. Quantitative assessment of these effects can be 
carried out using theoretical calculations of inclusive charged hadron production. 

4.4 Semi-inclusive recoil jet measurements 

Trigger hadrons lie within the charged-track acceptance |rj| < 0.9 and are selected in the intervals 8 < 
PT.trfg < 9 GeV/c, denoted by TT{8,9} and referred to as the Reference TT class; and 20 < PT.trig < 50 
GeV/c, denoted by TT{20,50} and referred to as the Signal TT class. 

Figure Q] shows semi-inclusive distributions (Eq. [3]) for recoil jets with R = 0.4, for the Signal and Refer¬ 
ence TT classes in pp collisions at ^/snn = 7 TeV and central Pb-Pb collisions measured by AFICE, and 
in pp collisions at ^Ann = 2.76 TeV simulated by PYTHIA. The distributions include all jet candidates 
in the recoil acceptance. 

Since p is the median energy density in an event, there must be jet candidates with energy density less 
than p, and which consequently have Pj j° t ’ ch < 0. The recoil jet distribution in the region pp j° t ’ ch < 0 
is seen to be largely uncorrelated with TT class in all cases, indicating that the yield in this region is 
dominated by combinatorial jets. In pp collisions the distribution in this region is narrow, indicating only 
small background density fluctuations. The predominant feature of the pp distributions is the strong de¬ 
pendence on TT class for pj j° ,ch > 0, with a harder recoil jet spectrum for higher ;?T,tng> as expected from 

the systematics of jet production in QCD. For Pb-Pb collisions the distribution in the region pj j°’ ch < 0 
is much broader, indicating significantly larger background density fluctuations than in pp collisions. For 
large and positive , the recoil jet distribution in Pb-Pb is strongly correlated with TT class, similar 

to pp collisions, showing that this region has significant contribution from the true coincident recoil jet 
yield. 

The integrals of the Pb-Pb distributions in Fig. Q] are 1.645 ±0.005(stat) for TT{8,9} and 1.647 ± 
0.009(stat) for TT{20,50}. This integral represents the average number of jet candidates per trigger 
hadron, both correlated and uncorrelated, and is seen to be consistent, within errors of a few per mil, for 
the two TT classes. Similar features have been observed in model calculations Il57l . Since central Pb-Pb 
events have high multiplicity, the recoil acceptance in each event is fully populated by jet candidates. 
Invariance of the integral with TT class therefore indicates that the number of jet candidates per trigger 
hadron is due largely to geometric factors, specifically the acceptance and jet resolution parameter R. 
This behavior is consistent with the robustness of the anti-kr algorithm against modification of jet struc¬ 
ture by soft particles from the underlying event l25l . Jet candidate distributions reconstructed using the 
anti-Lr algorithm for different trigger hadron kinematics appear to differ most significantly in the shape 
of the distribution as a function of p.j h - et , not in the total number of jet candidates per event. 

Based on these considerations we define a new observable, A reco ii, which suppresses the uncorrelated 
jet yield in a purely data-driven way. A reco ii is the difference between two semi-inclusive recoil jet 
distributions (Eq. [3]) for the Signal and Reference TT classes f57l . 
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p reco, ch ( GeV / C ) 




Fig. 1 : Semi-inclusive distributions of jets recoiling from a hadron trigger for two exclusive TT classes (Eq. 
0. for pp collisions at i/y = 7 TeV from ALICE data (top), pp collisions at i/s = 2.76 TeV from particle-level 
PYTHIA simulations (center), and central Pb-Pb collisions at yLviv = 2.76 TeV from ALICE data (bottom). All 
distributions are a function of /7^j°j Ch (Eq. 0. Distributions from data are not corrected for background fluctuations 
and instrumental effects. 
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The scale factor CR e f, which is within a few percent of unity, is discussed in Sect. 15.11 

The raw A reco ii distribution must be corrected for instrumental effects and for smearing of coincident 
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recoil jet energy by fluctuations of energy density in the underlying event. After corrections, A reco ii 
represents the change in the distribution of jets recoiling in coincidence with a trigger hadron, as the 
trigger hadron p-\ changes from the Reference to Signal TT interval. While this differential coincidence 
observable has not been reported previously, it is nevertheless well-defined in terms of perturbative QCD. 

We also extend Eq. @]to measure the angular distribution of recoil jet yield with respect to the axis defined 
by the trigger hadron momentum, in order to investigate medium-induced acoplanarity (3][28l (“inter-jet 
broadening”). The azimuthal correlation between the trigger hadron and coincident recoil charged jets is 
measured via 


d>(A<p) 


1 d 2 N jet 

Nttig d/?Tj et dA<p 


— CRef • 


/?T,trig£=TTsig 


1 d 2 N jet 

K^Pt^ A( P 


PT,trig€TTR c f 


(5) 


where the recoil acceptance for this observable is 7 t/2 < A (p < n. Normalization to unit i] is omitted 
from the notation for clarity. 

We quantify the rate of medium-induced large-angle scattering by measuring the integrated recoil yield 
at large angular deflection relative to A(p = n, defined as 

^ ^^lliresh 

£(A<Pthresh) = / dA<p[<E(A<p)], (6) 

Jn/2 

where the lower limit of the integration is set arbitrarily to 7 t/ 2. The upper limit excludes the main peak 
of the <E(A<p) distribution, \A(p — ji\< A<p t hresh> in order to measure the yield in the tail of the distribution. 
£ (Afresh) is measured as a function of Afresh. 

The distributions <f>(A<p) and £(A<p t hresh) likewise represent the change in the angular distribution of 
recoil jet yield, as the trigger hadron pj changes from the Reference to Signal TT interval. 

5 Raw distributions 

In order to ensure statistical independence of the recoil jet distributions for the Signal and Reference TT 
classes, each event is assigned randomly to one of the TT classes and is used only for its assigned TT 
class. The statistical reach of the analysis is optimized by assigning 80% of the events to the Signal TT 
subset and 20% to the Reference TT subset. This choice balances retention of the high-p^ j°’ ch component 
of the Signal recoil jet distribution with statistical precision of the Reference distribution in the region 
Pt j°t’ Ch < wid* the latter condition required to provide accurate normalization of the combinatorial 
background jet distribution. 

Events within each subset are then selected for further analysis if they contain at least one hadron within 
the pT.trig interval of their assigned TT class. If more than one hadron satisfying this criterion is found, 
one hadron is chosen randomly as the trigger hadron. For the Pb-Pb analysis there are 65k events with 
trigger hadrons in the Reference TT class and 22k in the Signal TT class. For the pp analysis there are 
74k events with trigger hadrons in the Reference TT class and 5k in the Signal TT class. 

5.1 Distributions of A reco ii 

Figure [2j left panels, show uncorrected trigger-normalized recoil jet distributions for R = 0.2, 0.4 and 
0.5 for both Signal and Reference TT classes. The right panels show the ratio of Signal and Reference 
distributions for each value of R. The error bars in Fig. [2] are statistical only, and are dominated by the 
statistics of the recoil jet yield in all cases. The statistical error due to trigger hadron yield is negligible. 
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Fig. 2: Uncorrected trigger-normalized recoil charged jet distributions for central Pb-Pb collisions, with Signal 
TT{20,50} and Reference TT{8,9}. Jets are reconstructed with the anti-Cr algorithm, constituent /?t, const > 0.15 
GeV /c, and R = 0.2, 0.4, and 0.5. Left: individual spectra. Right: their ratios. The red line shows a linear fit in the 
indicated region. Error bars show statistical errors only. Left-middle panel (R = 0.4) is identical to right panel in 

Fig.ffl 


Figure [3 left panels, also show the integrals of the distributions. The integrals are seen to be insensitive 
to TT class for a given R, with variations at the percent level or smaller, while the value of the integral 
depends strongly on R. These features are consistent with the geometric interpretation of the integral 
given in Sect. 14.41 


Table Q] shows the parameters resulting from the fit of a linear function to the ratios in the right panels 
of Fig. [2j in the region of pj j° ,ch where the distributions are largely uncorrelated with TT class. The 
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R 

PT,iet’ Ch m ran 8 e (GeV/c) 

Constant (cR e f in Eq. [4]) 

Slope (GeV/c) 1 

0.2 

[-12, 2] 

0.99 ±0.01 

0.004 ±0.002 

0.4 

[-20, 4] 

0.96 ±0.01 

0.004 ±0.001 

0.5 

[-26, 10] 

0.93 ±0.01 

0.002 ±0.001 


Table 1: Parameters from linear fits to ratios shown in right panels of Fig. [2] 


constant term of the fit, CR e f, is less than unity by a few percent, while the slopes exceed zero by about 2 
<7. The individual distributions vary by three orders of magnitude in this region. 

A value of CR e f below unity arises because the higher TT class has a larger rate of true coincident recoil 
jets, and the integrals of the distributions are largely uncorrelated with TT class. Larger yield at positive 
Pj j°t Ch consequently depletes the yield at negative and small positive values of j°' ch . 

Accurate subtraction of the uncorrelated component from the Signal TT distribution therefore requires 
scaling of the Reference TT distribution by CR e f, as indicated in Eq. [4]and Eq. [5] Scaling of the Reference 
TT distributions in the pp analysis by CR e f has negligible effect. 

Figure [3] shows A reC oii distributions for R = 0.2, 0.4 and 0.5. The left panels, which have linear vertical 
scale, show A reco ;i in the region of pj j°’ ch in which the scale factor CR e f is determined. A reco ;i is seen to 
be consistent with zero over the entire fitting range. These panels also show the close similarity of the 
shapes of the Signal and Reference distributions in this region. 

Figure Qj right panels, show A reco ii at positive Pj j° t ’ ch - where the Signal and Reference distributions 
diverge. This is the ensemble-averaged distribution of the trigger-correlated differential jet yield, but 
with measured pj j°’ ch not yet corrected for instrumental effects and fluctuations of the underlying event 
background. 

5.2 T>(A<p) and £ (A<p thresh ) 

The analysis of <f>(A<p) (Eq. [5]) and £ (A<p t h res h) (Eq. [6]) is the same as that for A reco ii in terms of event 
selection, track cuts, jet reconstruction, and jet candidate selection. For this analysis we only consider 
jets with R = 0.4 and 40 < pj j°’ ch < 60 GeV/c. 

Figure|4j left panel, shows the distributions of ( P(A(p) for TT{8,9} and TT{20,50} individually, and for 
TT{20,50}-TT{8,9}, illustrating the effect of the subtraction. 

Figure [4] right panel, shows the raw distribution of £(A<p t h re sh)> likewise for TT{8,9} and TT{20,50} 
individually, and for TT{20,50}-TT{8,9}. Since £(A<p t h re sh) is an integral over A(p beyond a specified 
threshold, care must be taken to ensure statistical independence of measurements for different values of 
the threshold. Each point in Fig. |4j right panel, is therefore generated from an exclusive subset of the 
data, with 10% of the data used for threshold values 0.1 and 0.2, 20% for 0.4, and 60% for 0.7. Subsets 
of unequal size are chosen to optimize the statistical errors. 

Due to the limited statistical precision of the data, correction of the raw distributions in Fig. [4] via unfold¬ 
ing for background fluctuations and instrumental effects is not possible. In order to compare the Pb-Pb 
distributions with a reference distribution for pp collisions, we therefore impose the effects of instru¬ 
mental response and Pb-Pb background fluctuations on the distribution calculated by PYTHIA for pp 
collisions at y/s - 2.76 TeV. The instrumental response, modeled by GEANT, is dominated by tracking 
efficiency and momentum resolution. The effects of background fluctuations are modeled by embedding 
detector-level PYTHIA events into real Pb-Pb events. Recoil jets are reconstructed from these hybrid 
events, using the same procedures as real data analysis. 
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Fig. 3: Distribution of raw A reco jj for A' = 0.2, 0.4, and 0.5, measured in central Pb-Pb collisions for Signal TT class 
{20,50} and Reference TT class {8,9}. Left panels: range of fit to extract CR e f, with linear vertical scale. 

Right panels: p'j j°’ ch range above fit region, with logarithmic vertical scale. Error bars show statistical errors only. 


12 














Semi-inclusive hadron-jet measurement in central Pb-Pb 


ALICE Collaboration 



s- 

<L 

w 


ALICE 

0-10% Pb-Pb ^ = 2.76 TeV 
AntL/tj charged jets, R = 0.4 


0 02 h 40 < p'; e “' ch < 60 GeV/c 

T,jet 


o.oi- 


f 


Hadron trigger 
• TT{20,50} - TT{8,9) 
O TT{20,50} 

D TT{8,9) 


■>- 


93 


0.1 0.2 0.3 0.4 0.5 


0.6 0.7 

A (p 


thresh 
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GeV /c. All error bars are statistical only. The vertical dashed line in the left panel indicates the acceptance limit 
for the A reC oii measurement. Points in the right panel are displaced horizontally for clarity. 


13 










Semi-inclusive hadron-jet measurement in central Pb-Pb 


ALICE Collaboration 


6 Corrections to A recoi i distributions 

Corrections to the raw A reco ii distributions for underlying event fluctuations and instrumental response are 
carried out using unfolding methods (581159 1. in which the true jet distribution T is determined from the 
measured distribution M using a response matrix. We denote by R t0 , the response matrix that incorporates 
all corrections, due to underlying event fluctuations and to instrumental response. R tot maps T to 

M(p %* et ), 


M(4jet) = *tot (p'Tjet)P t jet) >< T (Pj£ t)> ( 7 ) 

where p^jet' s particle-level charged-jet p-y and /Lj c j et is the detector-level or reconstructed jet p-\ . 

Precise inversion of Eq. [7]for non-singular R tot may result in a solution with large fluctuations in central 
values and large variance, arising from statistical noise in M(Pjj et ) EH- Inversion of Eq. [7] to obtain 
a physically interpretable solution is achieved via regularized unfolding, which imposes the additional 
constraint of smoothness on the solution. 

Input to the unfolding procedure uses jets in the range 20 < P^jet <100 GeV/c. The distribution in 
Eq. [4] provides a natural cutoff at low p^jet’ where the difference between central values of Signal and 
Reference distributions is smaller than the statistical error of the difference, so that imposition of a lower 
bound in this range is strictly speaking not required. However, in practice it was found that imposition 
of a lower bound at p^jet = 20 GeV/c, which is above the LO cutoff in terms of charged jet /;j! c j el , is 
needed for stable unfolding. This bound was kept as low as possible, to retain as much correlated signal 
as possible. The upper bound is set by the requirement that the highest p^jet bin has at least 10 counts. 

Correction for loss of jet yield in the excluded regions is earned out by applying a Pxjet-dependent 
efficiency £kj n , which is determined using PYTHIA simulations. £| {m is close to unity for all R in the 
analysis, over most of the range 20 < p^ t <100 GeV/c. Its value is £| an = 50% at p^ t = 20 GeV/c 
for all R, due primarily to detector efficiency, and £| (m = 70% at p^ t = 100 GeV/c for all R, due to the 
effects of momentum resolution and background fluctuations. The jet finding efficiency is 95% for p^t 
= 20 GeV/c and 100% for p^ t > 40 GeV/c, for all R. 

For the Pb-Pb analysis, the primary unfolding algorithm is an iterative procedure based on Bayes’ The¬ 
orem ia, as implemented in the RooUnfold software package l6H . Regularization is imposed by re¬ 
quiring only small variation between successive iterations, which occurs typically after three iterations. 
Closure of the unfolding procedure for A reco ji was tested in model studies with correlated spectrum and 
background fluctuations similar to those of this measurement 11571 . An alternative unfolding algorithm, 
regularized Singular Value Decomposition (SVD) (59], was used to estimate the systematic uncertainties. 

Both unfolding algorithms were also used for the pp analysis. In this case, the SVD algorithm was used 
to determine the central values, while Bayesian unfolding is used to estimate the systematics. This was 
found to be the optimal approach for the more limited statistics of the pp distributions. 

Both unfolding methods require initial specification of a prior distribution. For the Pb-Pb analysis, the 
prior is the A reco ji distribution for pp collisions at y/s = 2.76 TeV, calculated using PYTHIA (Perugia 10 
tune). For pp collisions at yfs = 7 TeV, the prior is calculated using PYTHIA (Perugia 10 tune (32l ). 

6.1 Correction for instrumental response 

The procedures to correct the jet energy for instrumental effects are the same as those described in (T71 . 
The dominant correction is due to tracking efficiency, with p-\ resolution generating the second-largest 
correction. 
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Corrections for instrumental effects are determined from simulations of pp collisions at \fs = 2.76 TeV 
generated by PYTHIA, together with detailed detector simulations generated using GEANT followed by 
event reconstruction. The lower tracking efficiency in central Pb-Pb collisions was modeled by randomly 
discarding additional detector-level tracks. The additional rejection factor was determined by comparing 
Hijing and Pythia efficiencies and corresponds to 2-3%, with weak px dependence | |2()1 . 

Jet reconstruction is carried out for each event at both the particle and detector level. The instrumental 
response matrix, Ra eU is generated by associating particle-level and detector-level jets whose centroids 
are close in (t],0), following the procedure described in i fTTI l. 

6.2 Correction for background fluctuations 

The adjustment of reconstructed jet energy by the estimated background density p • Aj et (Eq. [2]) accounts 
approximately for event-wise variation in the background level, which arises from variation in multiplic¬ 
ity within the 0-10% centrality percentile bin ll38l . The jet energy scale of the A reco ii distribution must 
still be corrected for energy smearing, due to local background energy density fluctuations relative to the 
median background density p. 

Background fluctuations Spj are measured by two techniques: the Random Cone method (RC) If38l . 
and a method in which model jets are embedded into real events f39l . The distribution of fluctuations in 
background energy for the RC method has RMS = 4.35 GeV /c for R = 0.2, 9.9 GeV /c for R = 0.4, and 
13 GeV/c for R = 0.5. The RC method is used for the central data points, with the embedding method 
used to assess the systematic uncertainty. 

The calculation of p (Eq. Q]) requires algorithmic choices that are not unique, notably the jet reconstruc¬ 
tion algorithm and the population of jets used for the median calculation. However, calculation of the 
response matrix for unfolding of background fluctuations incorporates the same set of choices. If all jet 
candidates are retained, without rejection based on p^j et , the effect of any systematic shift in JES due to 
p will be precisely counterbalanced by a shift of the same magnitude but opposite sign in the response 
matrix. This two-step JES correction, with event-by-event jet energy adjustment for event pedestal p -4 j et 
followed by ensemble-level unfolding of background fluctuations Spj, will consequently be independent 
of the specific algorithmic choices for determining p. 

In this analysis, the definition of p for the first step excludes the two hardest jet candidates from the 
median calculation (EqJTJ, while in the second step only jet candidates with p^jet > 20 GeV/c are used 
for unfolding. However, this rejection cut in the second step induces an implicit dependence on the 
specific definition of p. In order to assess this effect, the analysis was repeated with an alternative 
definition of p, in which all jet candidates are included in the median calculation in the first step. No 
significant differences were observed in the corrected recoil jet spectra. 

6.3 Cumulative response matrix 

For the Pb-Pb analysis, the cumulative response matrix R U)[ is the product of the response matrices for 
instrumental response and background fluctuations. To illustrate the magnitude of corrections to A reco ii 
from unfolding the raw distributions with R tot , we calculate the converse effect by convoluting the A reco ii 
distribution for pp collisions with R tot . Fig. [5] compares the particle-level A reco ii distribution calculated 
using PYTHIA for pp collisions at y/s = 2.76 TeV for R = 0.2 and 0.5 with their convolution with R^et and 
8p-\ separately, and R tot . The figure shows the ratio of the convolution over the unsmeared distribution. 
For R = 0.2, the effects of background fluctuations are small, and the net effect of R tot is due primarily to 
the instrumental response. For R = 0.5, the effects of background fluctuations and instrumental response 
offset each other to a large degree, with only a small net effect on the central values of the distribution. 
The distributions for R = 0.4 are similar to those for R = 0.5. Since the shape of the A reco ii distribution is 
similar in the pp and Pb-Pb analyses, the corrections in the two analyses will likewise be similar. 
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Fig. 5 : Ratio of A reco ;i distributions for pp collisions at yfs = 2.76 TeV, for R = 0.2 (left) and R = 0.5 (right). 
The numerator is convoluted with /?det and 8pj separately, and with R tol . The denominator is the unsmeared 
distribution. 


For the pp analysis, only the instrumental response was corrected using unfolding, i.e. R tot = P c iei • 

6.4 Other effects 

In this section we discuss other effects that do not warrant correction of the data. 

Since this analysis is based on a semi-inclusive observable, with normalization provided by the number 
of trigger hadrons measured offline, correction for online trigger efficiency (Sect. [2]) is not required. No 
significant difference in measured distributions was observed for events in the 0-8% and 9-10% centrality 
intervals. 

Tracking efficiency at high-pT is 80% (Sect. [2]), so that 20% of all trigger hadrons will not be ob¬ 
served. However, this tracking efficiency is uniform over the PT.trig range spanning both the Reference 
and Signal TT classes, so the loss in trigger statistics is unbiased in PT.trig- Since the measurements 
are trigger-normalized semi-inclusive distributions, the reduction in the observed trigger hadrons corre¬ 
sponds simply to a loss of events, and correction for this effect is not required. 

Section 14.31 presented considerations of trigger-jet energy loss in the interpretation of these measure¬ 
ments. A related but distinct effect is variation of Raa, the suppression of inclusive hadron yield in 
central Pb-Pb collisions, over the Pt, trig-interval of the Signal TT bin Ill3lll4ll . Such a variation can gen¬ 
erate different hard-process selection bias for the same hadron trigger cuts in pp and Pb-Pb collisions, 
even if trigger-jet energy loss effects in Pb-Pb are negligible. Using PYTHIA simulations, we estimate 
that this variation may generate an increase in A reco ji of 5% at PT,jet = 20 GeV/c and 15% at PTjet = 100 
GeV/c, but negligible change in the ratio of A reco ii in Pb-Pb with different R (see Fig.[[0]and discussion 
below). Such effects will however be included in theoretical calculations which incorporate quenching 
and accurately reproduce the measured pT-dependence of inclusive hadron Raa, and we do not correct 
the data for them. 

High-pf trig hadron triggers above a fixed p-\ threshold bias the event population due to correlation with 
the Event plane (EP) orientation, and bias towards more-central events. Both effects will bias the under¬ 
lying event density and its fluctuations in the recoil jet region. Note, however, that for the px.trig ranges 
of the Signal and Reference TT classes in this analysis, the second-order EP correlation amplitude V 2 ex¬ 
hibits no significant variation with hadron p-\ (v 2 ~ 0.01 for p-\ > 10 GeV/c 16211631 ). and the centrality 
bias is also invariant. The subtraction of the two distributions in A reco ii and <f>(A<p) thereby removes the 
effect of such background biases to a significant extent. Residual effects of these biases are assessed in 
Sect. |TJ and are included in the systematic uncertainties. 
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Multiple, incoherent partonic interactions (MPI) can generate both a trigger hadron and uncorrelated 
hard jets in the recoil acceptance, in the same Pb-Pb collision. A recent analysis of y-jet coincidences 
corrected for this background using a mixed event technique Il43l . Since the rate of uncorrelated hard 
interactions is by definition independent of pT.trig, the subtraction of the Reference from the Signal dis¬ 
tribution in Eq. [4] and Eq. [5] suppresses entirely the contribution of jet candidates from all uncorrelated 
sources, including jets found in the recoil acceptance arising from MPI. Correction for MPI effects is 
therefore not required, for all observables considered in this analysis. 

7 Systematic uncertainties 

7.1 Systematic uncertainties of A reco ii 

The systematic uncertainties for the distributions from Pb-Pb collisions are determined by varying pa¬ 
rameters and algorithmic choices in corrections for instrumental response and background fluctuations. 
For the A reco ji distribution for pp collisions at \fs = 1 TeV, systematic uncertainties are determined by 
varying the corrections for the instrumental response. 

The significant systematic uncertainties of the A reco ji distributions in Pb-Pb collisions are as follows: 


- Fit range for CR e f (Table [TJ: variation of limits for fit generates a variation in A reco ii of less than 

1 %; 

- Tracking efficiency: variation of A'j c t by changing the tracking efficiency by 5% generates a varia¬ 
tion in corrected A reco ii of 4% at pj h |( . L « 20 GeV/c and 15% at pj h |0t ~ 100 GeV/c; 

- Fragmentation model for instrumental response: determination of R & t using PYQUEN rather than 
PYTHIA. PYQUEN has large-angle radiation enabled and was tuned to LHC data lf64l . This gives 
a variation in A reco ii of 2% at pj h |eL pa 20 GeV/c and 13% at p^jet ~ 100 GeV/c; 

- Event plane and multiplicity bias: the trigger hadron yield and background fluctuation distributions 
are measured differentially in bins of azimuthal angle relative to the EP. The trigger hadron yield is 
found to be correlated with EP orientation, indicating non-zero elliptic flow. The response matrix 
is then obtained by weighting the azimuth-dependent background fluctuation distribution with the 
azimuth-dependent trigger hadron yield. The change in corrected jet yield with and without this 
weighting is less than 5%. The effects of the multiplicity bias are negligible; 

- Background fluctuations: using embedding rather than the RC method to measure background 
fluctuations (Sect. [672]) generates differences in the corrected A reco ;i distribution of less than 10%; 

- Variation in unfolding algorithm: termination of Bayesian unfolding after five rather than three 
iterations generates variations in A reco ii of~ 1% over most of the measured range. SVD unfolding 
yields A rec0 ii distributions that differ from the Bayesian-based corrected distributions by 1%; 

- Choice of unfolding prior: for Bayesian-based unfolding, the alternative priors are the A reco i] distri¬ 
bution for pp collisions at y/s = 2.76 TeV including a 10% or 20% relative energy shift, to model jet 
energy loss. For SVD unfolding, the alternative prior is the Bayesian-based unfolded distribution. 
The largest variation in the corrected A reco ii is less than 6% at all pj h j et ; 

- Spectrum binning and limits: variations of upper and lower spectrum limits generate variations in 
corrected A reco ,i of less than 3% at low pJjJj et , with negligible variation at high p^ h j et . Variation in 
choice of binning generates changes in corrected A reco ji that are less than 4%. 
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The products of weak decays make negligible contribution to /?T h j et because of the stringent track se¬ 
lection requirements of the analysis and the low material budget of the ITS and TPC. The systematic 
uncertainty of A reco ii due to the contribution of secondary vertex decays is less than 2%. 


Systematic uncertainty 

P?M = 25 GeV/c 

P C TM = 75 GeV/c 

Correlated 

Scale Factor CR e f 

(-1,40)% 

(-0,40.1)% 

Tracking efficiency 

(—6,+6)% 

(-16,+16)% 

Fragmentation model 

(—0,+2)% 

(-0,4-14)% 

EP bias 

(—2,+0)% 

(-3,40)% 

Uncorrelated or shape 

Background fluctuations 

(—0,+7)% 

(—0,+8)% 

Unfolding algorithm 

(-0,+l)% 

(-0,4-2)% 

Unfolding prior 

(—5,+0)% 

(—11,4-6)% 

Spectrum limits and binning 

(-1,40)% 

(-2,4-3)% 

Cumulative correlated uncertainty 

(-7,40)% 

(-17,+22)% 

Cumulative uncorrelated uncertainty 

(-5,4-7)% 

(-11,+11)% 


Table 2: Relative systematic uncertainties of A reco ii, for R = 0.4 jets and two values of Pj h j et for central Pb- 
Pb collisions at C^nn = 2.76 TeV, TT{20,50}-TT{8,9}. Uncertainties are expressed as negative and positive 
differences from the central values, with an entry of zero indicating negligible contribution. Uncertainties are 
classified as correlated and uncorrelated, as described in the text. 


Table [2]presents the significant systematic uncertainties for R = 0.4, at two values of p-y fcr Uncertainties 
are presented as the relative difference to the central values of the corrected A reco ii . We distinguish 
between correlated systematic uncertainties, arising from variations that generate a correlated change in 
the magnitude of the spectrum, and uncorrelated (or shape) uncertainties, arising from variations that 
preserve the integral but generate a change in the shape of the spectrum. Cumulative uncertainties are 
the quadratic sum of all correlated or uncorrelated uncertainties. Uncertainties for R = 0.2 and 0.5 are 
evaluated in a similar way. 

Similar systematic uncertainties were considered for the A reco ;i distribution in pp collisions at s/s = 7 
TeV as those discussed for Pb-Pb collisions. Uncertainty in tracking efficiency causes variation in A reco ii 
by 2% at PT.jet ~ 20 GeV/c and 9% at PT.jet ~ 60 GeV/c. The systematic uncertainty due to track 
momentum resolution is estimated to be 4% in the entire p-\ range. The shift of jet energy scale due 
to contamination by secondary particles and fake tracks causes a variation in A reco ii of less than 2%. 
Variations in the unfolding procedure, including change in the choice of unfolding algorithm, prior, and 
spectrum binning, result in A reco ji changes of ~ 5%. The cumulative systematic uncertainty is given by 
the quadratic sum of all individual uncertainties. 

7.2 Systematic uncertainties of 4>(A<p) and £ (A<p thl esh) 

The systematic uncertainties for the measurement of <4>(A<p) and £ (A<p t hresh) are presented in this section. 

- Since the yield of correlated hard jets decreases with increasing acoplanarity (i.e. increasing |A<p — 
n\), the scale factor CR e f in this region should approach unity. To assess this effect, CR e f is varied 
from its nominal values in Table Q] to unity. The resulting change in width of the uncorrected 
<l>(A<p) distribution (Eq© is 0.001 and the change in slope of the £(A<p t h res h) ratio (Fig. fl2l) is 
0.35, which are taken as the systematic uncertainties. 
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- Tracking efficiency less than unity will result in jets that are reconstructed from a subset of their 
charged track constituents, with consequent variation of jet centroid. However, the 5% relative 
uncertainty of tracking efficiency generates negligible variation in the width of the <f>(A (p) distri¬ 
bution and in the slope of the £(A<p t h re sh) ratio. 

- The EP bias due to the hadron trigger, discussed in Sect. 17.11 generates a change of 0.005 in the 
width of the <P(A(p) and 0.07 in the slope of the £ (Afpthresh) ratio. 


Systematic uncertainty 

Width of 4>(A<p) 

Slope of E (AtjPthresh) ratio 

Scale factor CR e f 

±0.001 

±0.35 

Tracking efficiency 

negligible 

negligible 

EP bias 

±0.005 

±0.07 

Cumulative uncertainty 

±0.005 

±0.36 


Table 3: Systematic uncertainties for the width of the <t>(A<p) distribution (Eqj8]i and the slope of the £ (A<p t h res h) 
ratio (Fig. [12] right panel). The cumulative uncertainty is the quadratic sum of all contributions. 

Table [3] shows all sources of systematic uncertainty for <4>(A<p) and £(A<p t h re sh), with the cumulative 
uncertainty given by their quadratic sum. Instrumental effects generate negligible uncertainty in the 
azimuthal correlations. 

8 Distributions for pp collisions at y/s= 2.76 TeV 

As noted above, comparison of Pb-Pb measurements to similar distributions in pp collisions at y/s = 
2.76 TeV requires calculations based on PYTHIA and NLO pQCD. In order to validate this approach, we 
compare PYTHIA and NLO pQCD-based calculations to ALICE measurements of A reco ji distributions 
in pp collisions at y/s = 7 TeV, using the data shown in Fig. Q] 

The NLO pQCD-based framework was developed initially to calculate the spin-dependent hadron-jet co¬ 
incidence cross section at y/s = 200 GeV ll27l . The calculation uses the DSS fragmentation function l65l 
and the CT10 NLO parton distribution function Il66l . We model hadronization by a shift in p-\ for parton- 
level jets ll67ll . with the magnitude of the shift determined by a fit to inclusive jet distributions lf68l . The 
resulting particle-level jet distribution is transformed to a charged-jet distribution by applying a response 
matrix calculated using PYTHIA. The systematic uncertainty of the resulting spectrum is estimated by 
independently varying the parton distribution function and the factorization and renormalization scales 
by a factor two. 

Figure |6j upper panels, show A reco ji distributions for R = 0.5, from ALICE data and calculations for pp 
collisions at y/s - 7 TeV. The lower panels show the ratios of these distributions to a function which 
parameterizes the ALICE data. The PYTHIA calculations for both tunes agree with the measurement 
within uncertainties; similar agreement is found for R = 0.2 and 0.4. The central values of the NLO 
calculation are above the measured data by about 20%, though the calculation is consistent with data 
within uncertainties for R = 0.5. The discrepancy in central values is larger for smaller R, reaching 
about 50% for R = 0.2, which is not consistent within systematic uncertainties. In pQCD calculations 
the difference between the parton and jet momenta involves an expansion in terms of log (A), whose 
contribution may be significant for small R f69l . Improved agreement between the NLO calculation and 
data for R = 0.2 may therefore be achievable using resummation techniques Il69l . 

Figure [7] shows the ratio of A reco ji distributions for R = 0.2 and 0.5 in pp collisions at y/s = 7 TeV. 
The measured ratio is compared with PYTHIA and NLO pQCD-based calculations. The grey boxes 
show the systematic uncertainty of the measured ratio, taking into account correlations of numerator and 
denominator. 
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Fig. 6: A reco ii distributions for R = 0.5, for ALICE data and calculations of pp collisions ^/s = 7 TeV. ALICE 
data, which are the same in both panels, are compared with calculations based on PYTHIA (left) and NLO pQCD 
(right). The green boxes in the right panel show the systematic uncertainty of the NLO calculation. The lower 
panels show the ratios of data and calculations to a smooth function fit to the data. 




Fig. 7: Ratio of A reco ;i distributions with R = 0.2 and 0.5, for ALICE data and calculations from pp collisions at 
yfs = 7 TeV. ALICE data, which are the same in both panels, are compared with calculations using PYTHIA (left) 
and NLO pQCD (right). 

The NLO calculation generates larger ratios than those observed in the data. A related observable, the 
ratio of inclusive jet production cross sections for R = 0.2 and 0.4 in pp collisions at y/H = 2.76 TeV, 
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has also been compared with pQCD calculations lf70l . This comparison shows that both hadronization 
corrections and perturbative effects that are effectively next-to-next-to-leading order (NNLO) in the in¬ 
dividual cross sections If68l are required for agreement. Perturbative QCD calculations to higher order 
than NLO are also needed to describe the ratio of A reco ii distributions presented here. 

In contrast, PYTHIA simulations agree within uncertainties with data, both for A rec oii at fixed R and the 
Arecoii ratio for two different values of R. These comparisons therefore favor PYTHIA calculations for the 
reference distributions at y/s = 2.76 TeV. PYTHIA combines a LO matrix element with a parton shower 
resummation of leading logarithmic terms of soft gluon radiation at all orders, leading to an improved 
description of data compared to a fixed order analytic calculation. 

9 Results 

9.1 A reco ii 



Fig. 8: A r ecoii for 0-10% central Pb-Pb collisions at yA^N = 2.76 TeV for anti-Cr jets with R = 0.2, 0.4 and 0.5. 
The vertical error bars are the square root of the diagonal elements of the unfolding covariance matrix, with the 
boxes indicating correlated and uncorrelated (shape) systematic uncertainties. 

Figure [8] shows corrected A reco ii distributions for central Pb-Pb collisions, for R = 0.2, 0.4 and 0.5. The 
shape of the distributions is approximately exponential, with larger per-trigger yield for R = 0.4 and 0.5 
than for R = 0.2. 

The R dependence of A reco ii is related to the distribution of jet energy transverse to the jet axis. Scattering 
of the parton shower within the hot QCD medium may broaden this distribution [|28]|7T]]. The magnitude 
of intra-jet broadening can be measured by comparing A reco ii distributions for Pb-Pb collisions with 
those for pp collisions, in which jets are generated in vacuum. We utilize two related observables for 
this purpose: (i) A/ aa , which is the ratio of A reco ii for Pb-Pb to that for pp collisions simulated using 
PYTHIA, for fixed R , and (ii) the ratio of A reco ii at two different R in Pb-Pb, compared with that in pp 
collisions. 

Figure [9] shows A/aa for R = 0.2, 0.4 and 0.5. Suppression of the yield of recoil jets in Pb-Pb collisions 
is observed, with similar magnitude for all R. 

Figure [I0]shows the ratio of A reco ii for R = 0.2 relative to A reco ii for R = 0.4 and 0.5, in central Pb-Pb and 
pp collisions. The systematic uncertainties of the Pb-Pb ratios take into account the correlated systematic 
uncertainties in numerator and denominator. The shape uncertainties are propagated independently in the 
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Fig. 9: A/\,\, the ratio of A reC oii in central Pb-Pb and pp collisions at yC = 2.76 TeV, for R = 0.2, 0.4 and 0.5. A reC oii 
for pp collisions are calculated using PYTHIA. 
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Fig. 10: Ratio of A reC oii for R = 0.2 relative to R = 0.4 (top) and to R = 0.5 (bottom), for central Pb-Pb (black) and 
pp collisions simulated using PYTHIA (red) at y/s = 2.76 TeV. 


ratio, since variation in their components induces different effects as a function of R. The distributions 
for Pb-Pb and pp collisions are seen to be similar, with no evidence for intra-jet broadening in central 
Pb-Pb collisions within the uncertainties. 

The CMS collaboration has reported a significant redistribution of energy within R < 0.3 for jets in 
central Pb-Pb collisions lf72l . potentially in contrast to Fig. [TO] However, that measurement and the 
one reported here cannot be compared directly. Modeling of the two measurements within the same 
theoretical framework is required for their comparison. 

Figures [9] and |T0] show that the recoil jet yield is suppressed, while the intra-jet energy profile is not 
changed significantly for/? < 0.5. We note in addition that the infrared cutoff for jet constituents (tracks) 
in this measurement is pr.const = 0.15 GeV/c, which strongly constrains the correlated energy within the 
jet cone that would not be detected by this measurement. 

Taken together, these observations are consistent with a picture in which there is significant in-medium 
transport of radiation to angles larger than 0.5 radians. This picture was initially suggested by a measure¬ 
ment showing that the energy imbalance of highly asymmetric jet pairs is compensated, on an ensemble- 
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averaged basis, by the energy carried by soft particles at large angles relative to the jet axis (23]. Also in 
this case, however, quantitative comparison of these measurements requires their calculation in a com¬ 
mon theoretical framework. 

The A re c 0 ii distributions in both pp and Pb-Pb collisions are well-described by an exponential distribution 
o< e~ p Tod/ 6 , with values of b around 16 GeV/c. Fig. [9] shows that A/aa has negligible dependence on 
PTjet * C)I ^ = 0-4 and 0.5 within 60 < ^Tjet < 100 GeV/c, which indicates that the values of b are similar 
within this PTjet range for the pp and Pb-Pb distributions. The value of A/aa in this region can therefore 
be expressed as the horizontal shift of an exponential distribution of fixed slope. For R = 0.5 in the range 
60 < Px h jet < 100 GcV /c, the suppression in A/aa corresponds to a shift in /;.j h |et of — 8 ± 2 (stat) GeV/c. 
In the scenario of negligible trigger-jet energy loss, this shift corresponds to the average partonic energy 
loss of the recoil jet population via energy transport to large angles, outside the jet cone. 

9.2 Azimuthal correlations 

Figure QT] shows the uncorrected <E>(A<p) distributions for central Pb-Pb data and pp simulations. As 
noted in Sect. 15.21 we compare the uncorrected <J>(A<p) distribution of Pb-Pb data to a reference distri¬ 
bution for pp collisions (PYTHIA, Perugia 2010 tune), modified by the background and instrumental 
effects expected for central Pb-Pb collisions. We recall that <F(A<p) suppresses the uncorrelated contri¬ 
bution from MPI, which otherwise would provide a significant background at large n — A (p. 



Fig. 11: <l>(A(p) distributions for 0-10% central Pb-Pb data (black circles) and pp collisions simulated by detector- 
level PYTHIA events embedded into central Pb-Pb events (red squares), at yATiv = 2.76 TeV. Jets have 40 < 
< 60, with p^ Ch not correc ted for background fluctuations and instrumental effects. The lines show the 
result of fitting Eq. [8] to the distributions, with the value of a from the fit as indicated. The error bars show 
statistical errors only. The Pb-Pb data points are the same as the solid circles shown in the left panel of Fig. [4] 

The absolute yield of the Pb-Pb distribution is seen to be smaller than that of the pp reference. This 
is consistent with the suppression observed for A/aa (Fig. |9j, which is the ratio of the integrals of the 
<f>(A<p) distributions over the range n — A<p < 0.6. 

The <t>(A(p ) distributions for Pb-Pb and pp collisions are characterized by fitting a function corresponding 
to an exponential plus a pedestal term If43l . 


/(A<p) = p 0 x e (A<p n)/a + pi, 


(8) 
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where the parameter a reflects the width of the distribution. The fit range is 2n/3 < A (p < n. The 
fitted values are Gpb-pb = 0.173±0.031(stat.) ±0.005(sys.) and (Tpythia — 0.164 ± 0.015(stat.), which 
are consistent within uncertainties. We find no evidence from this comparison for medium-induced 
acoplanarity of recoil jets with uncorrected energy in the range 40 < pj j°’ ch < 60 GeV/ c. 

The azimuthal distribution between a direct photon (pr.-y > 60 GeV/c) and a recoil jet ( pr jet > 30 
GeV/c) has been measured in central Pb-Pb collisions and compared to that from PYTHIA events 
embedded in a simulation of Pb-Pb collisions Il43l . Fits of an exponential function to these distribu¬ 
tions give similar values of a for central Pb-Pb and embedded PYTHIA, likewise indicating no evidence 
for medium-induced acoplanarity, though the values of a are larger than those for the analysis reported 
here. Comparison of the shape of the azimuthal distribution of di-jet pairs in central Pb-Pb data and 
embedded PYTHIA events has been reported j22ll23l . with indication of an enhancement in the tail of 
the distribution for central Pb-Pb collisions l22l . 
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Fig. 12: Left: E (A^) t h r esh) distributions for central Pb-Pb data (black circles) and the pp reference distribution 
(red squares), obtained by embedding detector-level PYTHIA events into real data. Right: ratio of the E(A<p t h res h) 
distribution between Pb-Pb data and the PYTHIA reference shown in the left along with a first-order polynomial 
fit. The error bars in both panels show statistical errors only. 


More detailed characterization of the change in the angular distribution with change in TT interval is pro¬ 
vided by £ (Atpthresh) (Eq. [6]), which measures the yield in the tail of the distribution beyond a threshold, 
nfl < A(p <71 — Afresh- Fig. Q21 left panel, shows the £ (A<p t hr e sh) distributions for central Pb-Pb data 
and the embedded PYTHIA reference. As discussed in Sect. 15.21 for the measurement of £ (A<p t h re sh) the 
dataset is divided into exclusive subsets, with each subset used for only one value of threshold, so that 
the data points in Fig. [12] are statistically uncorrelated. 

The relative contribution to £(A<p t h res h) of different physics processes may vary with A<p t | ires iv At suffi¬ 
ciently large angular deflection of the jet centroid, at a value of A<p t h re sh which has not yet been deter¬ 
mined, the yield is expected to arise predominantly from single hard (Moliere) scattering in the hot QCD 
medium |3]|28l. Figure [12] right panel, shows the ratio of the two distributions in the left panel. We 
utilize this ratio of absolutely normalized distributions in pp and Pb-Pb collisions to search for effects 
due to Moliere scattering. 

The value of the ratio at small A<p t hresh corresponds approximately to the A (p -integrated suppression in 
recoil yield in Fig. [9] while its dependence on Afresh provides a comparison of the shapes of the distribu¬ 
tions. This comparison is quantified by fitting a first-order polynomial function to the ratio of £ (A<p t hr es h) 
in the right panel. The fit gives a slope of —0.527 ±0.641(stat.) ±0.36(sys.), which is consistent with 
zero within uncertainties. If Moliere scattering were the only mechanism modifying the Pb-Pb distribu¬ 
tion relative to that of pp collisions the ratio at large A<p t hresh should be larger than unity; however, the 
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ratio is seen to be below unity at the largest measured A<p t hr e sh> indicating that other mechanisms have 
large effect in this region. We find no evidence from this measurement for medium-induced Moliere 
scattering. The uncertainty in this measurement is dominated by the statistical error, however, meaning 
that additional data, together with measurements at other jet energies and larger angular deviations, will 
provide more precise constraints on the rate of large-angle scattering in the hot QCD medium. 

10 Summary 

We have reported measurements of jet quenching in central Pb-Pb collisions at the LHC, using a new 
analysis method based on the semi-inclusive distribution of jets recoiling from a high -pj trigger hadron. 
Discrimination of coincident jet yield from background is carried out at the level of ensemble-averaged 
distributions with a trigger-difference technique, with no selection bias imposed on the recoil jet popula¬ 
tion. This approach enables measurement of the distribution of jets with large R and low infrared cutoff 
for jet constituents, over a broad range of jet energy. Distributions are reported for charged jets in the 
range 20 < p^jet < 100 GeV/c, for R = 0.2, 0.4 and 0.5. 

The differential recoil jet yield in central Pb-Pb collisions is suppressed relative to that in pp collisions 
by up to a factor two for 0.2 < R < 0.5. Together with the low infrared cutoff of this measurement, this 
indicates that medium-induced energy loss arises predominantly from radiation at angles larger than 0.5 
relative to the jet axis. The energy carried by this radiation, which is reflected in the magnitude of the 
spectrum shift under the assumption of negligible trigger-jet energy loss, is estimated to be 8 ± 2 GeV/c 
for charged jets with R = 0.5, in the range 60 < <jet < 100 GeV/c. 

The ratio of differential recoil jet yields with different R is similar for Pb-Pb and pp collisions. No 
significant medium-induced modification of the intra-jet energy distribution for angles R < 0.5 relative 
to the jet axis is thereby observed. 

The width of the azimuthal distribution of recoil jets relative to the trigger axis is measured to be similar 
in Pb-Pb and pp collisions for jets of 40 < Pj j° t ’ ch < 60 GeV/c. No significant medium-induced acopla- 
narity is therefore observed, consistent with findings from di-jet and direct photon-jet measurements. 
Large angular deflection of the recoil jet may be sensitive to the rate of single Moliere scattering in the 
hot QCD medium and provide a direct probe of its quasi-particle degrees of freedom. We observe no 
significant rate of such large-angle scatterings, though with limited statistical precision at present. These 
data, when combined with theoretical calculations, will provide guidance for the necessary precision to 
achieve discriminating measurements in the future. 
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